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Abstract 

DGLAP evolution equations may be presented in a form completely analogous 
to the Boltzmann equation. This provides a natural proof of the positivity of the 
spin-dependent parton distributions, provided the initial distributions at Qq are also 
positive. In addition, the evolution to Q 2 < Qq may violate positivity, providing 
therefore a 'time arrow'. The positivity condition is just | AP{j (z)\ < Pij(z) for z < 1 
for all types of partons, while the '+'prescription and terms containing 5(1 — z) do 
not affect positivity. This method allows one to complete immediately the existing 
proof of Soffer inequality at leading and next-to- leading order. 
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The positivity constraints play an essential role in the current analysis of various parton 
distributions, describing the nucleon spin structure. One should first mention here the 
Soffer inequality, putting an additional constraint for the quark transversity distribution 

i 

Ma;)|<g + (£)^[g(x)+Ag(x)]. (1) 

Also, the current next-to-leading (NLO) parametrizations [fj, |3| are chosen in such a 
way, that positivity for all helicity parton distributions is respected, i.e. \Af(z,Q 2 )\ < 

f(z,Q 2 ). 

One may wonder, to what extent the Q 2 evolution is compatible with the positivity 
constraints. In the present paper it is shown, that the relation is especially clear when 
the interpretation of the evolution equation as a kinetic equation is adopted. As an 
immediate consequence, the complete proof of Soffer inequality at leading order (LO) and 
next-to-leading order (NLO), is easily obtained. 

The probabilistic understanding of the evolution equations came already from the 
pioneering papers |4], |5|, 0. Actually, its standard form [] (for the time being, we shall 
confine ourselves to the non-singlet (NS) case) 

dq{x) _a ± f^Mpft, ( 2 ) 



dt 2it Jx y y' 
may be interpreted as a 'time' t = InQ 2 evolution of the 'particles' density q in the one 
dimensional space < x < 1 due to the flow from the right to the left, with the probability 
equal to the splitting kernel P. The key element of such an interpretation is the problem 
of the infrared (IR) singular terms in P, which was considered in detail some years ago ]7) 
(see also ||). The kinetic interpretation is preserved provided the '+' form of the kernel 
is presented in the following way 

P + (z) = P(z)-6(l-z) f 1 p(y)dy, (3) 

J o 

leading to the corresponding expression for the evolution equation 

This form has already a kinetic interpretation because the second term in the square 
brackets describes the flow of the partons at the point x @. It seems however instructive 
-'^For brevity, the argument t will not be written down explicitly in the parton densities. 
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H to make this resemblance even clearer by the simple change of variables z = y/x in 
the second term, which gives the equation the following symmetric form 

m = ^ [ r 1 dy M Pi i ) _r d ^Mp i v )] . (5) 

at 2,7i J x y y Jo x x 

It allows us to write it down in a way absolutely similar to the Boltzmann equation, 
namely 

dq(x) 

~HF = Jo dy ( a ( y ^ x "> q ^ ~~ a ( x ^ v)q( x ))i ( 6 ) 

where the one-dimensional analogue of the Boltzmann 'scattering probability' is defined 

as 

a s x 6(y> x) 

a(y ^ x) = —P(-) . (7) 

2n y y 

As a result, the known properties of the evolution equation receive the following simple 
interpretation: the cancellation of the IR divergencies between the contributions of the 
real and virtual gluons emission is coming from the equality of in- and out- flows when in 
both terms of @ one has y ~ x, following from the continuity condition for the 'particles' 
number. 

Also, the conservation of the vector current 

Jo dx ~dT~ = Jo dxdy ( a ( y ~^ x )q(v) ~ a ( x ~^ v)q( x )) = °> ( 8 ) 

comes from the integration of an antisymmetric function in a symmetric region. 

Morever, one can easily question the validity of the famous Boltzmann H-theorem. 
It happens not to be valid because of the violation of one of the key Boltzmann assump- 
tions, namely, 'molecular chaos' a(x — > y) = a(y — > x). The violation is due, in particular, 
to the (9— function, reflecting the causality properties, as it is well known. We shall also 
address this subject when discussing positivity, to which we are now ready to come. 

As the 'particles' density in the Boltzmann equation is positive by definition, this 
property is absolutely obvious. The negative second term in ([]) cannot change the sign 
of the distribution because it is 'diagonal' in x, which means that it is proportional to 
the function at the same point x, as in the l.h.s.. When the distribution gets too close to 
zero, its stops decreasing. This is true for both '+' and 5(1 — z) terms, for any value of 
their coefficient (if it is positive, it will reinforce the positivity of the distribution). 
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Let us consider now the spin-dependent case. For simplicity, we postpone the dis- 
cussion of quark-gluon mixing for a moment, but allow the spin-dependent and spin- 
independent kernels to be different, as they are at NLO. It is most convenient to write 
down the equations for definite parton helicities, which was actually the starting point in 
deriving the equations for the spin-dependent quantities || . Although the form, which we 
shall use, mixes the contributions of different helicities, it makes the positivity properties 
especially clear. So we have 

^P- = ^{P ++ i-) ® Q + (v) + p + -{-) ® ?-(y)), 

at Z7T y y 

^P- = ^(P + -(~) ® q + (v) + P ++ (~) ® <?-(*))■ (9) 
dt z7r y y 

Here P++(z) = (P(z) + AP(z))/2,P+_(z) = (P(z) - AP(z))/2 are the evolution 
kernels for definite helicities, and the shorthand notation for the convolution is adopted. 
As soon as x < y, the positivity of the initial distributions (q+(x, Qq), q-(x, Qq) — 0, or 
\Aq(x, Qq)\ < q(x, Qq)) is preserved, if both kernels P ++ , P+_ are positive, which is true, 
if 

\AP(z)\ < P(z), z<\. (10) 

The singular terms at z = 1 are not altering positivity, because they appear only in the 
diagonal (now in helicities) kernel P ++ (only forward scattering is IR dangerous). From 
the kinetic interpretation again the distributions q + , g_ stop decreasing, as soon as they 
are close to changing sign. 

Now to extend the proof to the quark gluon mixing is trivial. One should write down 
the expressions for the evolutions of quark and gluon distributions of each helicity 

rfg+M = a ( ) + pm^) 9 q _ {)) 

dt In y y 

+Pf + (-) ® G+{y) + Pf4~) ® G^y), 

y y 
^P = ^(p + -(-) ® ?+(*) + p ++ (~) ® 

dt In y y 

+Pf-{-) ® G + (v) + p£(~) ® G_(y), 

- ?(n G !(-) ® ? + (y) + ^(-) ® <?-(*) 



2vr ^y J * y 

+P^(-) ® G + (y) + P GG (-) ® G_(y)), 
1/ 1/ 



dG^x) = ^ {p G 1{ X ^ q+{y) + pG q{ X ^ { ) 

at In y y 

+P°E(*) ® G+ (y) + (-) ® G_(y)). (11) 

y 1/ 



If the inequality fllQ|) is valid for each type of partons [|l(| 



|AP^)| < « < 1; i,j = q,G, (12) 

all the kernels, appearing in the r.h.s. of such a system, are positive. Concerning the 
singular terms, they are again diagonal, now in parton type, and do not affect positivity. 
The validity of these equations at LO comes just from the way they were derived, as the 
(positive) helicity-dependent kernels were in fact first calculated in ref. ||. At NLO, the 



situation is more controversial [1C 



To conclude, the stability of positivity under Q 2 — evolution comes from two sources: 

i) the inequalities (0), leading to the increasing of distributions, 

ii) the kinetic interpretation of the decreasing terms. 

For the latter it is crucially important, that they are diagonal in x, in helicity and in 
parton type, which is related to their IR nature. 

It is interesting to note ||, that the 'back' evolution to Q 2 < (which is described 
by DGLAP equations as well, as soon as Q 2 >> Aq CD ), generally speaking, does not 
preserve positivity, as the regular and singular terms change their role. Namely, all the 
terms for z < 1, lead to the falloff of the distributions, while the singular terms make them 
increase. Starting from a 'pathological' positive distribution (say, increasing with x), it is 
possible to get the violation of positivity. This just means, that an arbitrary distribution 
cannot be obtained as a result of the evolution starting from lower Q 2 . As a result, the 
evolution kernels in fact impose some restrictions on the boundary conditions (mainly 
coming from the mentioned causality properties, defining some 'time arrow'), which are 
not easy to see in the original renormalization group equations for the moments. 

One should note in this connection, that the inequality analogous to fll2*|) for the 
moments of the splitting kernels (anomalous dimensions) is not sufficient for positivity, 
since the distribution with all the positive moments may be negative for small x. 

Let us now come to the evolution of Soffer inequality. According to the previous 
analysis it is straightforward to define the following 'super'-distributions 

Q+(x) = q+(x) + hx(x), 



Q-{x) = q + (x) - hi(x). 



(13) 



Due to Soffer inequality, both these distributions are positive at some point Qq, and the 
evolution equations for the NS case take the form 

dt 2n y y 

where the 'super'-kernels at LO are just 



(14) 
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C F Al + zf 



+ 35(l-z)\, 



2 L {i- Z y, 

p q_ [z) s ^g^-^r^)] (15) 

One can easily see, that the inequalities analogous to fll2|) are satisfied, so that both 
P+ + (z) and P+_(z) are positive for z < 1, while the singular term does appear only in the 
diagonal kernel. So, both requirements are valid and Soffer inequality is preserved under 
LO evolution. The extension to the singlet case is trivial, as the chiral-odd transversity 
distributions does not mix with gluons. Therefore, they affect only the evolution of quarks, 
and lead to the presence in the r.h.s. of the same extra terms as in ( p]) 

= ^(P? + Q ® Q + iv) + pU x ~) ® Q-(y) 

+Pl G (-)®G + (y) + Pf + (-)®G4y)), 

y y 
= ^(PU~) ® + ® Q-(v) 

dt Ltx y y 

+Pl G (~) ® G + (y) + Pg(-) ® GL(y)), (16) 

y 2/ 

which are all positive and free from singular terms, so positivity is preserved. 

Let us make a brief comparison with recent works in this direction. The LO analysis 
of JI1J emphasizes in fact the positivity of P+_(z) (whose validity is most nontrivial for 
valence quarks) and the role played by singlet terms, while the consideration of P+ + (z) 



and, in particular, of singular terms, are absent. The same remark applies to the recent 
NLO analysis |12| , where however, one can find the general reason for the absence of 
singular terms in non-diagonal kernel P+_(z) at NLO, which is of crucial importance for 
our proof. A partial NLO analysis was also performed in [I3|, where the specific role of 
the coefficient function was studied. The full analysis of the NLO case can now be treated 
in a very simple way. We just have to examine at NLO the 'super'-kernels considered 
above. This requires only the knowledge of P$ and Pjp (denoted AyP^ j n re f. [ 12 



and according to our previous analysis, Soffer inequality will be preserved at any Q 2 > Ql 
provided P2+\z) = [P$(z) + P { } {\z)}/2 and P^\z) = [PW(z)-Pf( 2 )]/2 are positive 
for z < 1, while the singular term is occuring only in the diagonal kernel. All these 
features appear clearly in the results of the numerical calculations shown in Fig. [I], where 
we have used the results of |T2[ which correspond to performing the NLO evolution in the 
MS scheme. 

One may also consider the analogue of the inequality ([J) for the complete NLO cal- 
culation including the coefficient function for the Drell-Yan (DY) process 0, 14]. This 
corresponds to the choice of the DY factorization scheme, when the DY cross-section is 
preserving its Born form ||14| , while the whole NLO result is assigned to the evolution of 
the parton distributions. 

Let us note, however, that Eq. ([I]) should not be necessarily valid for the whole NLO 
result. The reason is that it is not related to the positivity of the density matrix of a real 
particle, but rather to that of a quark in the hadron. This is in contrast to the well known 
positivity relation in polarized Deep Inelastic Scattering (DIS) |15| A 2 < y/R, where the 
observed virtual photon is considered instead. 

Although only the full NLO result is physically meaningful, this does not contradict 
the fact that there are relations, like ([!]), which are valid only partially, and should be 
considered like a restriction for the choice of the factorization scheme. 

For example, all the sum rules coming from the conservation of some operator (like 
fermion number or momentum conservation), generally speaking, are not valid for the full 
NLO result, and the factorization scheme should be chosen in order to preserve them. 

We show in Fig. |2| the calculations of the 'super'-kernels at NLO including the DY 
2 The correct NLO Q 2 evolution of hi, with no specific concern to the validity of Soffer inequality, has 
been also derived in lO. 
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coefficient functions from \12\ P^+^and V+- . One observes that while the non-diagonal 
one remains positive, this is not the case for the diagonal one v2+ which is negative in 
the small z region. However, if we now combine LO and NLO, and consider 

Vfl\z) = P$_(z) + ^T?l\z), (17) 

we see in Fig. |] that positivity is preserved even in the DY factorization scheme. 

In conclusion, we have presented the general proof of positivity of parton distribu- 
tions, based on the kinetic interpretation of the evolution equations. While the general 
reasoning is valid at any order of perturbation theory, the consideration of NLO requires 
to check the inequalities for the splitting kernels, for a particular factorization scheme 
|T0| , . Positivity is, generally speaking, not conserved for the 'back' evolution to lower 
Q 2 , providing a 'time arrow' and Soffer inequality is proven to be preserved by the LO 
and NLO evolutions. 
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Figure 1: The 'super'-kernels at NLO P++\z) and P+-\z) versus z. 



9 




10 




11 



